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Abstract 

Compared to purely atomic collisions, ultracold molecular collisions potentially support a much 
larger number of Fano-Feshbach resonances due to the enormous number of ro- vibrational states 
available. In fact, for alkali-metal dimers we find that the resulting density of resonances cannot 
be resolved at all, even on the sub-//K temperature scale of ultracold experiments. As a result, all 
observables become averaged over many resonances and can effectively be described by simpler, non- 
resonant scattering calculations. Two particular examples are discussed: non-chemically reactive 
RbCs and chemically reactive KRb. In the former case, the formation of a long-lived collision 
complex may lead to the ejection of molecules from a trap. In the latter case, chemical reactions 
broaden the resonances so much that they become unobservable. 
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I. INTRODUCTION 



The central conflict in ultracold molecular scattering is this: On the one hand, at ultralow 
temperature scattering observables are few in number, often limited to a single two-body loss 
rate, sometimes complimented by an elastic cross section, and limited to explicit information 
on only a small number of partial waves. On the other hand, the underlying dynamics that 
drives scattering consists of complex motion on a three- or four- (or more-) body potential 
energy surface (PES). This surface is moreover anisotropic, so that many more angular 
momentum states may contribute to scattering than the few represented by the asymptotic 
partial waves. How to properly distill the elaborate dynamics of the collision complex into 
observables remains an open question. Fj 

For ultracold scattering of alkali-metal atoms, the link between the PES and observables 
is cemented by the observation of Fano-Feshbach resonances. Here, the two-body PES's 
involved are comparatively simple, and the remaining undetermined parameters, consisting 
most simply of a pair of scattering lengths and a Cq coefficient, are used to fit data and 
to produce predictive models p]. Likewise, the observation of resonances should assist 
in interpreting collisions of molecules with light, isotropic partners such as helium [2], or 
perhaps even light molecules colliding with each other [3] . In such a case ab initio potentials 
are likely to represent something close to reality, and to be readily fine-tuned by fitting to 
resonances. 

In a recent paper we have begun to explore the role of resonant scattering on heavier 
species with highly anisotropic interactions, specifically, alkali-metal atoms colliding with 
alkali-metal dimers, at low temperature p£]. A main conclusion of Ref. [I] was that the 
density of states (DOS) of rotational and vibrational motion of the three-atom complex 
may be quite high, e.g., perhaps of order 1 per Gauss in Rb + KRb scattering. While 
these resonances may conceivably be resolvable experimentally, it is likely an impossible 
and unrewarding task to generate them explicitly from a PES. Rather, Ref. [I] adopted 
a statistical treatment of the resonances, asking what properties of the complex could be 
assessed on average. A unifying concept in this analysis was the mean decay width of the 
resonances, as given by the Rice-Ramsperger-Kassel-Marcus (RRKM) expression found in 
2 In cases where the molecules are chemically reactive at these temperatures, much more information could 

of course be extracted by state-selectively detecting the products of reaction, a task that has not yet been 

performed experimentally. 
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chemical transition state theory [5] 

Trrkm = (1) 

where p is the density of states in the vicinity of the collision energy, and N a is the number 
of open scattering channels. In ultracold collisions involving alkali-metal molecules, a large 
value of p and a small value of N Q (perhaps even N — 1) implies a dense forest of very 
narrow resonances. 

In the present paper we extend this analysis to collisions of pairs of alkali-metal dimers. A 
main finding is that the DOS for the four-atom complex is vastly larger than for the 3-atom 
complex, so that the resonances so formed cannot be resolved at all, even on the T = 0.1 — 1 
pK temperature scale of experiments. In this "iiberresonant" regime, all observables become 
averaged over many resonances, effectively bypassing the inherent intricacy of the complex. 
The resulting non-resonant cross sections are then in principle actually easier to compute 
and interpret than the atom- molecule case. 

We apply this idea to two cases. One case is RbCs, which is not chemially reactive at 
ultralow temperature, and for which therefore iV = 1 in its absolute, that is, ro- vibrational 
and spin ground state. In this case, the resulting extremely narrow resonances imply long 
complex lifetimes, potentially on the order of experimental times. This means not only that 
some fraction of the molecules remain "invisible," hidden inside four-body complexes, but 
also that the complexes, upon colliding with another molecule, can be ejected from the trap, 
leading to an unwelcome delayed-three-body loss mechanism. In this article we provide 
estimates of the loss rates implied by this mechanism, including the effect of electric fields. 
A second example is afforded by KRb molecules, which remain chemically reactive even at 
ultracold temperature [Hj. In this case N Q includes all possible channels of the products of 
reaction, and is quite large. Thus the resonance width implied by ([I]) is far larger than the 
mean resonance spacing, and resonances are expected to be unobservable in the loss rates. 

II. THEORETICAL MODEL 

We consider collisions of diatomic molecules AB (where A and B denote alkali-metal 
atoms) in their 1 H electronic ground state, their v = vibrational ground state, their n = 
rotational ground state, and some nuclear spin states /4M4, I B M B , assumed to be decoupled 
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in a magnetic field. We pay attention to the nuclear spins in order to completely specify the 
state, and to properly account for Bose/Fermi symmetrization, but they play little other role 
in the theory we describe below. Moreover, let \LMl) denote the partial wave of the incident 
channel, describing the relative orbital angular momentum of the molecules. An important 
quantity in the theory is then the total angular momentum (exclusive of the nuclear spin) 
J = ni + 112 + L. Since we consider only asymptotic states with n\ = n 2 = 0, the value of 
J is identical to the partial wave L in a given incident collision channel. 
Introducing the shorthand notation 

|a> = S\ l E,v = 0,n = 0,lAM A ,I B M B )i 
x\ 1 Z,v = 0,n = 0, I A M A , I B M B ) 2 

x \LMl) (2) 

(where S denotes the appropriate symmetrization for bosons or fermions) , the collision cross 
sections can be written in terms of the scattering matrix elements S a > a , 

0"a^a' = ^2 l 1 — ^a'o| 2 ^- (3) 

LM L L' M' L 

k is the wave number of the colliding molecules and A accounts for their indistinguishability, 
that is, A = 2 if they are in identical states and A = 1 otherwise. The indices a, a' summarize 
the quantum numbers of AB in the incident channel, and are extended to include the product 
channels in the case of reactive collisions. Even in an electric field, the projection of the 
total angular momentum onto the field axis is conserved. 

Following Ref [3], we construct a scattering theory that incorporates both a high den- 
sity of resonant states of the collision complex, and threshold effects relevant to ultralow 
energies. This is achieved by combining multichannel quantum defect theory (MQDT) with 
the methods of random matrix theory. In doing so, we exploit the conceptual difference 
between the spin channels \a) that describe physics at large interparticle separation R; and 
the numerous resonant states of the complex, denoted that differ by rotational and 
vibrational quantum numbers from a. The key feature of MQDT is that one only needs to 
provide the reactance matrix K^ r which is defined at a "matching radius" R m that defines 
the boundary between short- and long-range physics. The MQDT formalism as outlined in 
Refs. [31 IB] accounts exactly for the wave functions for R > R m and directly yields the 
physical scattering matrix g phys via standard algebraic procedures. 
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As in our previous work, for R > R m we assume simplified long-range interactions of the 
form 

where E a is the threshold of the ath channel, which may depend on a magnetic field B. Here, 
m r is the reduced mass of the scattering partners and C@ is their van der Waals coefficient, 
which is taken to be isotropic in this model. These potentials are used to calculate the 
relevant MQDT parameters from which the cross sections are ultimately constructed. We 
will see below how to account for nonzero electric fields. 

The short-range i^-matrix is constructed according to the dictates of random matrix 
theory jH] as 

It is indexed by the N a asymptotic channels a, but is influenced by the myriad (i.e., N 3> N a ) 
resonant states /i. The input parameters for the resonant scattering theory, Eq. ([5]), are the 
zero-order positions of the resonances and the coupling elements W afi to the asymptotic 
channels. Within our statistical framework, and W afJ- are taken as random variables 
based on the Gaussian Orthogonal Ensemble (GOE) jH [9j. By employing such a model, we 
assume that the collision complex corresponds classically to a long, chaotic trajectory that 
ergodically explores a large portion of the allowed phase space. 

The GOE is in turn specified by the mean resonance width. It was determined in JJ] that 

a reasonable approximation for this width is the RRKM result itself, 

N a 

r = - — = Trrkm, (6) 
2np 

where N a is the total number of asymptotic channels in the v = 0, n = ground state 
manifold. Further narrowing of the resonances due to the Wigner threshold laws is accounted 
for within the MQDT theory. 

Thus the resonance model is completely specified by the density of states p. We estimate 
the DOS in the same way as in Ref. j^. Namely, we posit a set of approximate potential 
curves 

h 2 L (L + 1) 

Vlj(R) H ^ ° R2 h E rv (v cl , n c i, v c2 , n c2 ). (7) 

Here Vlj(R) is a Lennard- Jones potential with the correct C§ for the molecule-molecule 
interaction, and tuned to a depth equal to the binding of the A 2 B 2 complex relative to the 



AB + AB threshold. This potential is augmented by a partial wave L c of the complex, and 
by a threshold energy E TV corresponding to ro- vibrational excited states of the molecules in 
the complex. Key to our DOS approximation is that all possible states that preserve the 
total angular momentum J = n c i + n C 2 + L c and conserve energy are included. Although 
the total J is limited to a few values as dictated by the incident partial wave L of scattering, 
the angular momentum quantum numbers of the complex can span into the hundreds (see 
below). 

Having identified all such relevant potentials ([7]), we compute their bound states lying near 
the incident threshold, and by counting them determine the DOS. The complete DOS thus 
constructed assumes ergodicity, i.e., that all states not forbidden by conservation laws are 
actually potentially populated. However, this assumption can be adjusted by, say, reducing 
the maximum value of orbital angular momentum L c used in the estimate. 

III. FOUR-BODY DENSITY OF STATES 

We estimate the DOS as described above and in Ref. [I], for two prototypical ultracold 
molecules: RbCs [TT] and KRb [121 [13] • To construct the Lennard- Jones potentials in Eq. (J7|) 
for these species, we use the C 6 coefficients from Ref. |14j], and potential depths of 800 cm -1 
for (RbCs)2 [2] and 2779.6 cm -1 for (KRb)2 [15] . To compute the ro- vibrational spectrum 
E rv we employ the empirical potential of Ref. [16] for RbCs, and that of Ref. [17] for KRb. 

The resulting ro- vibrational DOS for several total angular momenta J is reported in Table 
|TJ In Ref. [I] the possibility for processes that change the nuclear spin were considered, but 
we do not do so here; thus the table counts only the ro-vibrational density of states. Also 
shown is the mean lifetime of the collision complex, estimated as r = 2irhp. These estimates 
assume that all states of allowed angular momentum defining the complex, L c , n cl , n c2 in the 
vicinity of threshold can actually be populated. Even relaxing this assumption and reducing 
L c , the DOS remains quite high, as seen in Fig. [I] 

An exemplary elastic cross section for RbCs molecular collisions in the absolute ground 
state (v = n = 0, Mq s = 7/2, M^b = 3/2 for both molecules) is presented in Fig. |2} Shown 
are cross sections for s-wave [yellow (light gray)] and rf-wave [orange (gray)] scattering. This 
figure covers an energy range of twice the van der Waals energy -E v dw — ft 3 (2m r )~ 3 / 2 C 6 ] 
and contains 10,000 energy points, evenly spaced on a logarithmic grid. At this resolution, 
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TABLE I. Ro- vibrational densities of states for M = and corresponding one-open-channel RRKM 
lifetime ro = 2irhp; for N > 1 one has r = tq/N q . Regarding the magnetic dipole moment, in the 
case of (fermionic) KRb we assumed for odd J that the both molecules are in their lowest state 
(Mk = —4, M R b = 3/2); for even J one molecule is considered to be in this lowest state, the other 
in the next higher one (-Mk = —3, M R b = 3/2). For (bosonic) RbCs it is the other way round: for 
even J they are in the same state (MRb = 3/2, Mq s = 5/2) and for odd J in the next higher one 
(M Rb = 3/2, M Cs = 7/2). Values used for the magnetic moments are [TO]: ( 40 K 87 Rb) M K = -4, 
M Rb = 3/2, Mmag = 2.84kHz/G; M K = -3, M Rb = 3/2, Mmag = 3.08kHz/G; ( 87 Rb 133 Cs) M Rb = 
3/2, Af Cs = 5/2, /i mag = 3.50kHz/G; M Rb = 3/2, M Cs = 7/2, ^ mag = 4.07kHz/G. 
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most of the s-wave resonances are resolved, and the cross section frequently approaches the 
unitarity limit, 8n/k 2 . For rf-wave scattering, these resonances are not resolved so well. It is 
clear from the figure that at typical ultracold temperatures T = 0.1 — 1 K, these resonances 
can never be resolved. 

IV. INFLUENCE ON SCATTERING OF NON-REACTIVE MOLECULES 

The extremely high density of states estimated in the previous section implies a striking 
feature of molecule-molecule cold collisions. Namely, molecules that meet on resonance may 
become lost in the complex for times on the order of many milliseconds, comparable to the 
time scales of a typical experiment. In this section we formulate a set of rate equations 
accounting for this occurrence, using RbCs as an example. The rate equations describe 
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maximum angular momentum 

FIG. 1. (Color online) Ro- vibrational density of states of RbCs + RbCs collisions as a function 
of the maximally allowed end-over-end angular momentum L c of the collision complex, for a total 
angular momentum J = — 4, M = (solid, long-dashed, dashed, dashed-dotted, and dotted line, 
respectively). The scale on the right axis gives the corresponding RRKM lifetime of the collision 
complex according to tq = 2irhp. 




collision energy (K) 

FIG. 2. (Color online) Elastic s-wave [yellow (light gray)] and d-wave [orange (gray)] cross section 
of RbCs + RbCs. Incident channel is the absolute ground state (v = n = 0, Mq & = 7/2, Mr^ = 3/2 
for both molecules). 

three separate events: i) a pair of RbCs molecules meet and stick together, thus temporarily 
transforming into four-body complexes, with number density n c ; ii) The complexes decay 
back into molecules on a time scale set by the mean lifetime of the resonant states; and 
iii) during the lifetime of the complex, another RbCs molecule can collide with it, leading 
almost certainly to trap loss. We deal with each of the three parts of this process in the 
following. 
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A. Molecule-sticking rate 



For RbCs in its absolute ground state, the number of open channels is exactly N — 1. 
This circumstance automatically places resonant scattering in the limit where, on average, 
resonance widths are smaller than the mean resonance spacing, and resonances do not over- 
lap. Thus only some fraction of the collision events lead to long-lived resonances, albeit very 
long-lived ones. We model the sticking process by ascribing to it a cross section which is 
zero away from resonance, but which contains resonances at the appropriate DOS and width 
distribution. Such a cross section is in fact afforded by the elastic cross section a el (E), which 
is easily computed from our statistical MQDT formalism. 

The rate at which the complex-forming collisions happen is given by a thermally averaged 
rate constant, here distinguished by the partial wave L, M L considered: 

PCX) 

K^\T)= / a?, ML (E)vf(v)dv, (8) 
Jo 

where "mm" stands for "molecule-molecule," and 

™- (9) 

is the Maxwell-Boltzmann distribution for the relative velocity for a given temperature T of 
the initial molecular sample. 

When the mean resonance spacing is far less than the temperature, as we assume, then 
these many resonances are averaged over. We can therefore replace the strongly-varying 
cross section by its mean value, taken over each of the isolated resonances separately, and 
averaged over the mean spacing d = p^ 1 between resonances: 

i pE+d/2 

<m l {E) = -J <M L {e)de. (10) 

a J E-d/2 

This amounts to saying that only a fraction of collision energies, approximately T(k)/d, are 
on resonance and can lead to large sticking times, where T(k) is the mean resonance width 
in the vicinity of energy E = h 2 k 2 /2m r . Note that the resonant cross sections scale as the 
unitarity limit, oc l/k 2 , whereas resonance widths T(k) oc k 2L+1 , leading to a threshold law 
Kmm L ^ 2L f° r ^e sticking rate. 

More quantitatively, we make use of the simple algebraic structure of the MQDT formal- 
ism, in the ultracold limit and for a single channel (N a = N a = 1) the elastic cross section 
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reads 

rn-A 4?r A{E ' L) ' fin 

l,Ml \ i lA 2m r E[K^{E)]-^ + A{E,Lf' 



where K sv = —^^2 IJi W^/(E — E^). In the ultracold limit, the energy dependent MQDT 
parameter A(E,L) can be written down explicitly [8], 

^■ L)1/2 - r(f + !) r(, + |) ^(7l|l) ' (12) 



-Rvdw = (2m r CQ/h 2 ) 1 ^ i is the van der Waals length scale. In deriving Eq. (11) we employed 
the ultracold limits Q — > and 77 — > of the remaining MQDT parameters [8]. In the 
vicinity of a resonance at E and replacing the short- to long-range couplings by their 



average, ttW 2 —T/2 jl], Eq. (11) becomes 



Assuming additionally that A(E, L) is approximately constant within the range of a single 
resonance, Eqs. ( flQp3| ) yield the mean cross section at collision energy E, 



A-7T 



We therefore identify the rate constant for collisional sticking as 



poo 

K^ Ml) (T)= / at ML (E)vf(v)dv (15) 
Jo 



A r(f + i) 2 r(L + l)^ ( ^ T) • 



Interestingly, this expression agrees exactly with the inelastic rate constant derived for 
scattering in the presence of rapid loss due to chemical reactions [18j, modeled by assuming 
unit loss probability at each collision energy. The effect of averaging over a very large 
number of very narrow resonances has produced a cross section that is equivalent to full 
absorption at every collision energy, modified by the appropriate threshold laws. This is a 
tremendous simplification: rather than even attempt to deal explicitly with real potential 
energy surfaces and the many resonances they engender, we are able to cut immediately to 
the observable consequence, namely, temperature-dependent sticking probabilities. 

Armed with this insight, we can immediately extend the model to nonzero electric fields, 
assuming that the field significantly affects only long-range physics, and molecules reaching 
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FIG. 3. Quenching rate constants of two indistinguishable bosonic polar 87 Rb 133 Cs molecules as 
a function of the induced dipole moment for a temperature of T = 1 ^uK. The solid line represents 
the Ml = contribution, the dashed line the Ml = 2 (equal to the Ml = —2) contribution. 

small R vanish with unit probability. This problem can be solved exactly as in Refs. [TBI HH] ■ 
The resulting rate constants for our example of RbCs collisions are reproduced in Fig. [3} 
References [HJ [19] predict that rate constants for loss in partial wave L scales as d iL+2 for 
induced dipole moment d. Hence, for small dipole moments s-wave scattering prevails. In 
Fig. [3] the rate constant in the upper curve shows an initial rise oc d 2 for s-waves, until it 
saturates at around d ~ 0.25D |20j. There is then a second rise, owing to the rapid increase 
of loss rate in the o?-wave channel, which dominates the loss beyond d ~ 0.6D. For scattering 
with orbital angular momentum component Ml = 2 (lower curve), the ci-wave rise is still 
apparent, but there is, of course, no s-wave contribution at smaller dipole moment. 




B. Mean lifetime of the complex 

The resonant complexes formed in molecule-molecule collisions will eventually decay back 
into pairs of molecules. The lifetime of the complex at a given collision energy can be 
quantified by means of the time delay [2TH23] . 

dS 

r delay = 2h— . (17) 

Here 5 is the eigenphase sum, that is, the sum of the inverse tangents of the eigenvalues of 
the i^-matrix. Employing the same approximations as in deriving Eqs. 
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13114), the time 



delay close to a resonance at E reads 



A(E,L)T/2 

Tdelay ~ (E-E y + [A(E,L)T/2}^ (18) 



and therefore the mean time delay becomes 



rE +l/(2p) ^ 

tdelay = P / 2H 7^ dE = 2?r V ( 19 ) 

JE -l/(2p) d& 



This is just the lifetime of the resonant complex as predicted by the RRKM theory, Eq. (|6j). 
Just as we need not consider individual resonances in the high-density limit, neither do we 
need to consider their individual lifetimes - another simplification. We therefore define, for 
each partial wave L, a decay rate of the complexes, jl = ^deiay' which follows immediately 
from the DOS in Table HI 



C. Rate equations 

We are now in a position to formulate the rate equations for the ultracold gas of RbCs 
molecules. Denote by n m the number density of these molecules and by n Cj L the num- 
ber density of the transient four-body complexes formed from initial partial wave L of the 
molecule-molecule scattering. Because the molecule-molecule scattering rates and 
the decay rates are different for different L, we explicitly add together the different con- 
tributions, as if they were independent. We assume that molecule- complex collisions are 
s-wave dominated and field independent, and hence described by a universal rate of the 



form Eq. (16), with L = and appropriate values for the reduced mass and Cq. Rate 
equations that describe the sticking of two molecules to form the complex, the subsequent 
decay of the complex, and demolition of a complex due to collision with another molecule, 
are given by 

> Ml) + 2 lL n c , L - K mc n m n c , L I , (20) 



n, 



L \ M L 



n c ,L = \n 2 m ^2 K mr^ L) ~ lLn c ,L - K mc n m n c>L . (21) 

M L 

Here Kmm lL ^ and K mc are the molecule-molecule and molecule-complex collision rate con- 
stants. As shown in Table [TJ different total angular momenta J lead to different densities of 
states and therefore different lifetimes of the collision complex (recall that J is identical with 
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the incident partial wave L of the collision) . This is accounted for in Eq. ( 20 ) by allowing the 
formation of different, independent collision complexes with densities n C) L, each possessing 
its own decay rate 72,. The complexes are populated according to the molecular collision 
rate Kmm 1 ^ for the given partial wave L as extracted from Fig. [ij The molecule-complex 
collision rate constants are considered equal for all complexes. Moreover, we assume that 
different Ml give rise to the same DOS and therefore to the same lifetime 7^. 

The time-dependent molecular density n m (t) resulting from a numerical integration of 



Eqs. ( 20pT ), starting from an initial molecular density n , is presented in Fig. |3](a). Results 



are shown for two different electric field strengths (0 and 10 kV/cm, respectively). Initially, 



when n c {t) <C n m (i), Eq. (20) is dominated by the loss due to complex formation at a 
rate — K mm nl, where K mm = Ml Kmva 1 ^ is the total molecular loss rate. After some 
time if this initial, fast decay turns over into a slow decay due to lossy molecule-complex 
collisions. Some insights can be gained by setting K mc = for the moment, that is, no lossy 
molecule-complex collisions. The resulting molecular density n^? (t) is shown in Fig. 4^a) as 
dotted (assuming only s-wave collisions) and dashed (c?-wave collisions) lines. In this case, 
the solution rim\t) reaches an equilibrium, 



ns?w - Vi+ y-r'" 1 . < 22 > 

'-''mm/ 

- -*L (23) 

as t ^> t' . The timescale t', on which the initial linear decay turns to reach this dynamical 
equilibrium, can be extracted from the analytic solution (t) as 

t' = . 2 (24) 
A/7 2 + K miI an Q 

This time is indicated in Fig. |4j^a) by arrows, assuming purely s-wave collisions for zero field 



and ci-wave collisions for lOkV/cm. Moreover, by inserting Eq. (23) into Eq. (20) one finds 
an expression for the slow final decay, 

n m (t) w [1 + aK mm K mcl - l n\(t - tj] ~ 1/2 , (25) 

where n m (tx) = n\ for some time t\ > t' at which the long-times behavior has already been 



reached, a acts as an empirical correction factor that accounts for Eq. (23) not reaching the 
dynamical equilibrium quite yet. 
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FIG. 4. (a) Time evolution of the RbCs molecular density n m {t) for zero field (upper curves) and 
an applied electric field of lOkV/cm (lower curves); s- and d-w&ve wave collisions are considered. 
Dotted (s-wave) and dashed (d-wave) lines are computed without particle loss due to complex- 
molecule collisions, that is, K mc = in Eqs. ( 20|21 ). The initial molecular density is no = 10 10 cm~ 3 



at a assumed sample temperature of 1 fiK. The arrows indicate the timescale t' as given by Eq. (24). 
(b) Evolution of the ratio of molecules bound in complexes to free molecules. The same set of 
parameters as in panel (a) is used. 



The time evolution of the molecular density is vastly influenced by external electric fields, 
scaling as d iL+2 for dipole moment d and partial wave L [18]. As a result, for our example 
in Fig. |4](a), the molecular density after its initial, fast decay is almost cut in half for fields 
> lOkV/cm. Even the field-free case is not free of losses due to complex formation; however, 
over 90% of the initial density is retained. Hence, in spite of not being chemically reactive, 
ultracold RbCs may still manifest substantial loss, which will become faster in an electric 
field. 

This is emphasized by Fig. |4](b), where we show the ratio of total molecules bound in 
complexes to free molecules, for the same set of field strengths as in panel (a). As expected 
from Fig. Qa), for zero electric field only a very small fraction of the molecules is bound in 
complexes. This changes drastically, however, once the field is turned on. For lOkV/cm, 
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half the molecules are trapped inside collision complexes at any given time. 



The apparent loss of molecules depends on the magnitude of any applied electric field, 



but also on the initial density n of molecules, cf. Eq. (22). In Fig. |5|we show the results for 
uq = 10 11 cm~ 3 , that is, increased by one order of magnitude compared to Fig. El Now, even 



zero electric field leads to significant molecule loss. For lOkV/cm, after only half a second 
fewer than 20% of the initial molecular density is retained. 



Experimental data on loss as in Figs, gjg can strongly constrain the parameters of 
the model At short times, the initial decay is given by — i^ mm ^o- With the initial density 
uq usually well known, one can therefore extract K mm straightforwardly. At intermediate 



times, n m (t) turns from its initial drop to its long-time behavior. This timescale, Eq. (24), is 
proportional to the lifetime of the complex. Hence, from experimental data, one could infer 
at least an estimate of the complex' lifetime. Finally, at long times the decay of the molecular 



density is well fitted by Eq. (25), from which in turn the molecule-complex collision rate 
constant K mr can be extracted. 
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V. INFLUENCE ON SCATTERING OF REACTIVE MOLECULES 



The situation is completely different for a species that is chemically reactive at zero 
temperature, such as KRb. In this case, transition state theory dictates that the number 
of open channels N Q includes also the product channels. These appear to be shockingly 
numerous, considering that the exoergicity of ground state KRb collisions is only 10.4 cm -1 
[6]. Even within this small energy release, it is possible, in principle, to produce K2 molecules 
with rotation quantum number up to = 13, or Rb2 molecules up to nR& 2 = 20, or any 
energetically allowed combination. Moreover, the products can have any reasonable partial 
wave angular momentum L K2Rb2 of the products about each other, provided that the total 
angular momentum J = n K2 + n Rb2 + L K2Rb2 is conserved. 

Again assuming ergodicity in all degrees of freedom, we obtain N Q by simply counting 
all possible exit channels consistent with conservation of angular momentum and energy, 
constructing molecular levels from the potentials in [T71 EH EH]. The result is a vast num- 
ber of possible channels, which grows rapidly as a function of total angular momentum J. 
Accounting for all these possibilities, the resulting number of possible product channels are 
listed in Table [TT] for various total J but fixed M = 0. No enhancement due to nuclear spins 
is considered here. 

Also shown in the table is the corresponding RRKM decay rate into product channels. 
This width far exceeds the mean level spacing (by a factor of N a , in fact), and renders the 
individual resonances unobservable. In fact, in this limit one expects collision cross sections 
to exhibit Ericson fluctuations, on a scale comparable to T itself. Inasmuch this width is 
already of order /zK (or ~ 10 Gauss in magnetic field), it is unlikely that any structure will 
be seen at all. Again, we are back to the simpler situation of studying non-resonant cross 
sections. 

Indeed, the occurrence of many exit channels implies that the decay rate of the complex 
is extremely rapid, so rapid that the states of the complex are not significantly populated 
at all. An alternative way to view this circumstance is to note that in the statistical theory 
the probability of chemical reaction is N /(N + 1) ph 1, whereas the probability of elastic 
scattering back to the single initial channel is 1/(Nq + 1) < 1. Thus the scattering leads to 
unit probability of chemical reaction, as posited in Refs. [T8ti20| [261429] . In fact, Ref. [28] 
provides a universal analytic expression of the inelastic rate constant for p-wave scattering 
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TABLE II. Number of KRb + KRb -> K 2 + Rb 2 product channels along with the predicted RRKM 
width, Trrkm = N /2-Kp, of the resonances. In calculating the latter, we use the densities of states 
provided in Table [TJ 



J number of channels Irrkm^K) r RRK M(G) 

45055 2.21 7.78 

1 131239 2.15 7.27 

2 213521 2.11 7.42 

3 291901 2.06 6.97 



with unit reaction loss probability, for identical fermionic molecules, 

K^f = Un^axikafA. (26) 

Here, a = 2 7 r J R vdW /r(l/4) 2 and a x = a r(l/4) 6 /[12vrr(3/4)] 2 are length scales of s- and 
p-wave scattering from a pure C§ potential. 

In Fig. [6] we compare representative inelastic p-wave scattering cross sections to the unit 



loss prediction Eq. (26). Rather than immediately employ the full number of open channels 
N Q for this case (which is technically challenging even within our simplified model), we 
emphasize the trend for ever-larger N a . Thus the solid black line shows an exemplary cross 
section for only N Q = 100 open channels. In this case Ericson fluctuations occur on a sub-yuK 
scale and are seen in the spectrum. However, even increasing the number of open channels 
to N Q = 1000 (red line) almost completely washes out these fluctuations. Moreover, this 



result shows almost perfect agreement with the simple model in Eqn. (26). We find similar 
good agreement for different realizations of the statistical spectrum. 

We conclude from this result that the realistic N a , which is much larger still, will certainly 



lead to a featureless loss spectrum given by Eqn. (26 ). Thus the statistical model as deployed 
here vindicates the models in Refs. [1814201 [26TT291. 



VI. SUMMARY 



For collisions of alkali-metal dimer molecules, we have found that the density of ro- 
vibrational states is enormous, far too large to probe individual resonances even at the 
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FIG. 6. (Color online) Inelastic p-wave scattering cross section mimicking ultracold KRb collisions 
with N Q = 100 (light gray) and N Q = 1000 [red (dark gray)] product channels. The dashed line 



represents the analytic prediction for unit loss probability, cf. Eq. (26). 



sub-/iKelvin energy resolution afforded by ultracold temperatures. Because of this circum- 
stance, resonant collision rates are always averaged over many resonances, and the theoretical 
description of scattering is greatly simplified. Thus broad general conclusions can be drawn. 
For the case of chemically reactive molecules, the formation of a resonant state ensures that 
the atoms have ample opportunity to find their way into the product channels, at least for 
reactions that are assumed to be barrierless. This in turn leads to essentially unit probability 
of reaction in each collision event, consistent with interpretations of recent experiments in 
ultracold KRb gases. 

Strikingly, even molecules that are not chemically reactive at zero temperature, in the 
presence of this vast number of resonances, behave as if they were chemically reactive, at 
least transiently. These molecules are capable of sticking together for a finite lifetime, which 
is dependent on the density of states. The longer this lifetime is, the more likely that the 
molecules bound in resonant complexes will be struck by other molecules and lost. Contrary 
to expectation, it may therefore be necessary to shield even non-reactive molecules from 
collisions by confining them to ID lattices and immersing them in electric fields [29H31] . 

A main quantitative uncertainty in the results described here is whether the full density 
of ro- vibrational states is actually populated in a collision, which may not be the case [35J. 
If it is not, then the time during which the molecules are stuck together reduces, and the 
loss rates may not be as great. Thus measurements of loss may provide direct insight into 
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the four-body dynamics of the molecule-molecule complex. 
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